A binuclear atom — a special type of close bound state 
between proton and heavy atom 
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It is established that a bound state of a proton with a heavy Thomas — Fermi atom should exist. 
On the one hand, the electrons of the atom screen the proton field. This decreases the repulsion 
force between the proton and the nucleus. On the other hand, the attraction force between the 
proton and the electrons is directed towards the gradient of the electron density, i. e. towards the 
nucleus. For instance, for Z = 80 both forces become equal at approximately 0.6 a where a is the 
Bohr radius. The corresponding minimum of the proton potential energy is in the region of negative 
energies (attraction) that can be of the order of several tens of eV. We propose to call such a system 
a binuclear atom. 

In contrast to the molecules where a coupling with a hydrogen atom is due to an essential modi- 
fication of one or several states of the outer electrons the formation of a binuclear atom is a result 
of collective response of the whole system of inner electrons to the screened potential of a proton 
that is well inside the electron system of the heavy atom. The variation of the wave function of each 
electron can be considered as a small perturbation. The bound state is formed as a result of joint 
action of a large number (of the order of Z) of perturbed inner electrons. The important problem 
concerning the accuracy of our calculation is discussed. 

PACS numbers: 31.15.Ne, 31.15.Gy, 31.15.Md, 31.15.Bs 

INTRODUCTION 

The purpose of the present paper is to demonstrate that a proton can be bound to a heavy atom with the charge 
of the nucleus Z 3> 1. The principal idea of the paper can be formulated as follows: when a proton approaches the 
nucleus its potential is screened by the atomic electrons. 

As is shown by Teller [l[ (see also Ref. Q), two Thomas — Fermi (TF) atoms cannot form a bound state. In the 
present paper we consider what can be looked upon as an opposite case. The principal difference of this case and the 
problem considered by Teller can be formulated as follows. If one considers two heavy atoms the Teller theorem states 
that such a system has a lower energy than the same system where the atoms are merged into a molecule. However, 
in the present paper we consider a situation where one of the interacting items is a heavy atom [Z ^ 1) that can be 
described by the TF theory whereas the second item is a proton with Z = \. Our purpose is to prove that in this 
case (not covered by the Teller's theory) a bound state can be formed. The screened potential of a proton can be 
considered as a perturbation for the electrons of the TF atom. The bound state is formed as a result of joint action 
on the proton of a large number of the perturbed inner electrons of the atom. The vibrational energy of the proton, 
i. e. the energy of its vibrations near the equilibrium position is much smaller than the binding energy. Therefore 
one can consider the proton in the S-state to be positioned on a sphere of a fixed radius. 

There is a number of papers where bound states between a positron and an atom are considered, see for instance, 
Refs. 0,0, [H| and the references therein. The ability of positrons to bind to a number of atoms is now well established. 
The attractive electron-positron interaction leads to formation of a cluster that includes a positron in the outer valence 
region of the atom. 

In contrast to the problem of a bound state of a positron the methods and results of these papers cannot be directly 
applied to the situation considered here because of a great difference between the positron and proton masses. In the 
binuclear atom treated in the present paper the distance between the heavy nucleus and the proton is smaller than 
the Bohr radius a. It means that in this case there is a quite different physical situation where the proton interacts 
practically entirely with the inner electrons of the atom. 
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PRELIMINARY CONSIDERATIONS 

To visualize the physics of this phenomenon we start with analysis of a subsidiary problem. We will consider a 
Thomas — Fermi screening of a proton field by a homogeneous electron system having average density N neutralized 
by immobile positive charges. We have 
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5N=—i^e^, V2,^ + qV = 0, </.(i?) = -cxp(-gi?), (1) 



p is the Fermi momentum, q is the reciprocal screening radius of a proton y, §40. Now, 
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7vi/3 = i:^±7vi/3. (2) 
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Here e is the elementary charge, m is the electron mass, N is their concentration, 5N is its variation due to the 
electrostatic potential if. The energy of electrostatic interaction of the screened proton potential situated at i? = 
with the electron density N is 
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E = -e^N / (fR- exp(-gi?) = -j:-ae^N^/^. (3) 

J R 3^/3 



Let us assume that is a slow (as compared with 1/q) function of R: N = 7V(R). Then the electron system will 
act on the proton with the force F given by 



dE /7r\4/3 
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Its direction is determined by VrA^ 



F = -^ = 2g) ^ ae'N-'^HR)VnN. (4) 



CALCULATION OF THE INTERACTION POTENTIAL 



Now we will treat the problem we are interested in, namely a proton in the field of electrons and a heavy nucleus. The 
electron system will be treated within the model of Thomas — Fermi atom (see, for instance, Ref. §70; fl|). This 
model was proposed long ago but quite recently an interest to its application was revived (for instance Refs. [gUlolfllt). 
According to the model half of the electron charge of the atom is within a sphere of the radius 1.33a.Z~^/3 where a is 
the Bohr radius. The electron density enhances as the distance from the nucleus becomes smaller. Let the distance 
between the nucleus and the proton be of the order of a. Its potential should be screened by the electrons. As a 
result, the repulsion between the proton and nucleus becomes much smaller. As the electron density increases while 
R decreases the electrons should attract the proton. The force of attraction cannot be calculated using Eq. ^ as 
VrTV now is not small. The calculations below show that for Z = 80 the force of attraction by the electrons and 
the force of repulsion by the nucleus are counterbalanced for R « 0.6 a. This is an indication for the possibility of a 
bound state formation. We will be interested in the bound state with the lowest energy, i. e. the lowest S-state. 

We wish to emphasize that there is a principal difference between the unique state treated in the present paper 
and an ordinary molecule of a chemical compound containing a hydrogen atom. In such compounds the coupling of 
a hydrogen atom is due to an essential modification of one or a few outer electron states whereas here we consider a 
special effect that is a collective change of the whole distribution of the atom's inner electrons by the proton potential. 
In ordinary molecules a coupling with a hydrogen atom is due to an essential modification of one or several states of 
the outer electrons. In the present paper we consider a special effect that is a collective response of the whole system 
of inner electrons to the screened potential of a proton. In regard to each electron the action of the potential can be 
considered as a small perturbation. The bound state is formed because these perturbations are summed up. Thus 
the perturbation theory is an adequate approach to this problem. Usually the perturbation theory cannot describe a 
bound state formation. However, the present case is unique as the bound state is formed as a result of united action 
of a large number of perturbed electrons. 

Consider the screening of a proton by the electrons of a heavy atom. We assume the distance between the proton 
and nucleus to be ^ a. In other words, the proton is well inside the atom. This is why we propose to call such a 
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system a binuclear atom. We will make estimates of the energy of such system within the TF atom model. In the 
first order of the perturbation theory the variation of electron density SN is 

<5iV(r) = 4 V V(A:|^|z) 11™(0,x)W(0,x)V'n;WV'n'i'M. (5) 

% k 

Here i = nlm are the occupied states, k = n'l'm' , n' and/or I ^ l'\ 9, x are the spherical angles, i/jni is the radial 
quasiclassical (WKB) electron wave function in the potential <I>(r) of the TF atom. Here and henceforth we use the 
atomic units e ^ h = m = I. 

Let us transform Eq. ([5]) into such a form that it gives a local relation between SN and ip. 6N being given, ip is 
determined by the Poisson equation. Inserting the relation between SN and ip we get a selfconsistent equation for ip. 
The variation of the electron energy in the first order of the perturbation theory is 

Ssi = {i\(p\i). (6) 

The perturbation for every electron is the proton potential plus the potential produced by SN, i. e. the screened 
proton potential. In ^ and © , there is a contribution of a selfacting electron. It gives a contribution to the screened 
potential (p. This contribution is however negligibly small as the electron number Z is large. The total energy variation 
of the system AEi is the sum of expressions ([5]) over all the particles including the nucleus 



AEi = Zip{0) - J N{r)ip{r)(fr. (7) 
Let z-axis join the nucleus and the proton. Then ip{r) does not depend on the azimuthal angle and can be expanded 



as 



ip{r) = ^ ips{r)Ps{cosl 



(8) 



where Ps is a Legendre polynomial. As N{r) has a spherical symmetry, only (poif) contributes to Eq. ([7|). Now we 
will expand SN like Eq. ([5]) and integrate Eq. ([5]) over 9 using the well-known identity 



\Y,^\' = {4n)-\2l + l) 
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As a result, Eq. (fS)) takes the form 

SNoir) = TT-i ^(2/ + 1) E E ('''\^o\n 
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The radial quasiclassical function of the state (n, /) with the energy near the Fermi level is 



; ^ \ -1 -1/2 

Wni(r) = anir Pi cos 



(pi + SniPi ^)dr' + Ci 



, = 2<i>Kr) ^ 2$(r) - - I ; + i 



,1/2 



(9) 
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(11) 



{Ci are constants of the order of 1). In Eq. ([TTjl the approximation {jpf + 2e„()"'' ~ pi + £niP^'^, is used. It is valid 
for \eni\ <C where i"/ are the characteristic values of ^i{r), a^i ~ (2/7r)(ie„;/dn — ^Tj^^, Ti is the period of the 
classical motion at the Fermi level, ri (as well as r2) are the classical turning points. 
Introducing the variable 



we expand ip{){r) into Fourier series 



(12) 
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Discarding the quickly oscillating part of the product of functions given by Eq. (fTT|) one gets 



5No{r) = V ^^^^ V ^ou 



de' J de 



cos[(e' - s)ti] 
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The integral over e' can be presented as a sum of integrals over the intervals [— , 0] and [0, oo] . The contribution of 
the first integral vanishes as the change of the variables e^e' changes the sign of this expression. The conditions e < 
and e' > are fulfilled for the first (5- function and e e [—2TruT^^,0]. In the expansion Eq. only the harmonics 
where <C are physically important. Integrating over e' with regard of (5-function and over e in the interval 

2'Ki>Ti~^ ,0, and taking into account Eq. p2| one gets 



SNoir) 



2Tr^r'' 
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(14) 



r is in the classical region for the quantum numbers I satisfying the condition <&i(r) > 0. Integration over I from up 
to the maximal value, corresponding to <E';(r) > gives 



SNoir) = TT-Mr)Mr), P = (2$)^/^ 



(15) 



Eq. ((T5)) coincides with the first equation ([T]) . Thus the relation between SN and (p is preserved also for nonhomo- 
geneous electron density. Eq. ([T5|) gives the Poisson equation 
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r dr'' 
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with the boundary conditions 



r = i?- 
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(16) 
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The jump of the derivative is due to the presence of the proton at the distance R from the nucleus. For small 
values of r we solve these equations numerically. For r > R one can use the eikonal approximation that is valid at 
d^/dr <C 2p^q. Using the TF variable x ~ ^R, 7 = ^^/■^/0.885 one can see that these solutions (but not the last 
boundary condition at x = 7i?) are independent of Z. For x < jR xipo{x) — CP{x), P{x) is given in Table 1. 

Table 1 
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P{x) 
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P{x) 


x 


P{x) 


0.2 


0.207 


1.4 


2.25 


2.6 


6.71 


0.4 


0.437 


1.6 


2.79 


2.8 


7.78 


0.6 


0.703 


1.8 


3.40 


3.0 


8.97 


0.8 


1.01 


2.0 


4.09 


3.2 


10.3 


1.0 


1.37 


2.2 


4.87 


3.4 


11.7 


1.2 


1.78 


2.4 


5.74 


3.6 


13.3 



For a; — > 0, P = X. For x > jR 

xipo{x) = bq~^^^ exp 



- J q^dx 1 , ql = 1.5- x-^^l'^X^I'^, $ = ZXjj 



(18) 



The constants c and h are determined by continuity of function x^po and the jump of derivative of function pT|) at 
X — ^R. 

Making use of the well known relation for the Thomas — Fermi atom between N{r) and ^(r) and the values of $ 
given in the tables one can calculate A£'i Eq.([7]) for various R. Note that an additional correction is needed (that, 
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FIG. 1: Fig. 1:1 — the energy of the system proton — heavy atom {Z=80) as a function of the distance R between the proton 
and the nucleus. 2 — the potential For _R > 1 one can expect that curve 1 will tend to $(-R). A possible transition 

curve is indicated by the broken line. 



however, will be small as compared with AEi). We have considered the screening in a neutral atom. The total 
electron charge surplus near the proton due to the screening is —1. Indeed, the total potential of the proton and SN 
Eq. tends to zero at r — * cx3 faster than r^^. The corresponding correction AE2 should not exceed the energy 
of one-electron ionization of the atom. In the TF atom model this energy appears to be lower as compared to the 
experimental values. Such a discrepancy is due to the fact that in this model such properties are essentially related to 
the distances r 3> 1 where the TF approximation fails. The average (over a number of different atoms) experimental 
value of the ionization energy is about 0.3 atomic units. The variation of the energy AE — AEi + AE2 is presented 
in Fig. 1. Here we have assumed AE2 = 0.3 a. u. 

Let us analyze the approximations used to obtain AE. First, this is the one-electron approximation. The char- 
acteristic value of the kinetic energy of the electrons should be bigger than the potential energy of their Coulomb 
interaction, i. e. '3> N^^"^. For the atoms where Z ^ 1 this condition is satisfied. Second, this is the quasiclassical 
approximation. It is satisfied in the interval < r < 1. The solution (|18p decays exponentially for r > R. This is 
why the integral ([7]) is determined by the interval r < _R. 

We have calculated the energy AE for Z=8Q and R < 1. For R > 1 our theory ceases to be applicable. [When 
calculating the integral ([7]) the interval of integration has been taken from to 1 where the integrand is well defined] . 
N{r) decreases with r; <po(^) increases for r < R and decreases for r > R. The integrand varies slowly for r < R and 
rapidly decays for r > R. Thus for R < 0.8 the discarded additional energy is small. 

Fig. 1 shows that the energy of the system proton — heavy atom {Z « 80) with regard of screening appears to be 
negative at > 0.4 a. u. and has a minimum at i? w 0.6 a. u.(« 0.3 A). One can expect that at the periphery of the 
atom AE can become positive (potential barrier). 

The distance between the bound proton and the nucleus (w 0.3 A for Z w 80) is much smaller than the characteristic 
radius of the valence electron wave functions. This means that the chemical behavior of such a system after an extra 
electron has joined it should be similar to that of an atom of the number Z+1. 
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CONCLUSION 

The experimental search of the predicted phenomenon should be probably performed with plasmas or using a source 
of low-energy protons. The energy release in the course of proton trapping might be accompanied by Auger-processes 
and X-ray emission. The probability of a proton capture might be bigger for a heavy atom in a condensed matter 
(i. e. for a target) than for a vapor. The described method is applicable for estimates of the binding energy of two 
protons with a heavy atom. If the average distance between the protons is of the order of 0.3 A the Coulomb energy 
of their repulsion is smaller that the coupling energy obtained in the present paper. Our estimates indicate that a 
possibility exists for trapping of a-particles and "'"-mesons. 

The TF model describes the average atomic characteristics. It cannot describe their individual properties, such as 
their periodicity. This point concerns also our estimates. 

A few words about the accuracy of our calculation. In general the accuracy one can expect of the TF atom model 
is about 10 — 15 per cent. Calculating the potential curve one deals with a difference of two large quantities. This 
means that the difference of these quantities should be bigger than 15 per cent of each of them. Near the minimum 
of the potential curve these two quantities differ by 30 per cent. This means that the accuracy of our calculation 
is sufficient. However, it is desirable to apply in future more advanced methods to obtain a more detailed value for 
the potential minimum position as well as the form of potential curve near the minimum. This should facilitate the 
experimental search of the binuclear atoms. 
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Figure caption 

Fig. 1. 1 — the energy of the system proton — heavy TF atom {Z=80) as a function of the distance R between 
the proton and the nucleus. 2 — the potential ^{R)- For i? > 1 one can expect that curve 1 will tend to $(i?). A 
possible transition curve is indicated by the broken line. 



